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Abstract 

Motivated by microscopic traffic modeling, we analyze dynamical systems which have a piecewise 
linear concave dynamics not necessarily monotonic. We introduce a deterministic Petri net extension 
where edges may have negative weights. The dynamics of these Petri nets are well-defined and may 
be described by a generalized matrix with a submatrix in the standard algebra with possibly negative 
entries, and another submatrix in the minplus algebra. When the dynamics is additively homogeneous, 
a generalized additive eigenvalue may be introduced, and the ergodic theory may be used to define 
a growth rate under additional technical assumptions. In the traffic example of two roads with one 
junction, we compute explicitly the eigenvalue and we show, by numerical simulations, that these two 
quantities (the additive eigenvalue and the growth rate) are not equal, but are close to each other. 
With this result, we are able to extend the well-studied notion of fundamental traffic diagram (the 
average flow as a function of the car density on a road) to the case of two roads with one junction 
and give a very simple analytic approximation of this diagram where four phases appear with clear 
traffic interpretations. Simulations show that the fundamental diagram shape obtained is also valid for 
systems with many junctions. To simulate these systems, we have to compute their dynamics, which 
are not quite simple. For building them in a modular way, we introduce generalized parallel, series and 
feedback compositions of piecewise linear concave dynamics. 

I. Introduction 

The main purpose of this paper is to explore some new points on dynamical systems that have 
emerged with the study of microscopic traffic modeling using minplus algebra and Petri nets. 
These points are the following : 
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• Standard Petri nets or linear maxplus algebra have not enough modeling power to describe 
the dynamics of a traffic system as simple as two roads with a crossing and with the simplest 
vehicle move. 

• A remedy to this difficulty is to introduce possibly negative weights on Petri net arcs (not to 
be confused with negative places or negative tokens as in [41]) or to use nonlinear minplus 
dynamics that are not mono tonic. 

• A consequence is that all known standard methodology based mainly on dynamic program- 
ming is ineffective. 

• Nevertheless, it is sometimes possible to use the standard ergodic theory, and/or to compute 
a generalized eigenvalue for the subclass of additively homogeneous systems and get nice 
qualitative results on the system. 

• When it is not possible to derive analytical results, we can do numerical simulations. In 
general, a microscopic modeling of a realistic system has a very large number of state 
variables difficult to manage. However, for the construction of the model, a generalized 
matrix calculus using minplus algebra can be developed for the class of system described 
here (having piecewise affine concave dynamics). 

All these points are illustrated on the traffic application for which we have obtained a new result: 
a good analytic approximation of the macroscopic law called the fundamental diagram in the 
case of roads with junctions. 

The traffic on a road has been studied from different points of view on a macroscopic level. 
For example: 

• The Lighthill-Whitham-Richards Model [34], which is the most standard view, expresses 
the mass conservation of cars seen as a fluid: 



where ip(x,t) denotes the flow at time t and position x on the road; p(x,t) denotes the 
density; and / is a given function called the fundamental traffic diagram. For traffic, this 
diagram plays the role of the perfect gas law for the fluid dynamics. The diagram has been 
estimated using experimental data, and its behavior is quite different from standard gas at 
high density. We obtain an idea of the diagram's shape in the subsection titled "traffic on 
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a circular road." 

Daganzo in [13] using the variable N(t,x) counting the cumulated number of vehicles 
having reached the point x at time t, and remarking that tp = d t N and p = —d x N has 
interpreted the equation ip = f(p) as an Hamilton- Jacobi equation when / is concave. 
The other equation d t p + d x (p = telling that d tx N = d xt N. Then, in [14], Daganzo and 
Geroliminis generalized this to a network by postulating the existence of an Hamiltonian 
in this case also and by trying to approximate it as an infimum of linear constraints that 
he obtains by physical considerations on the street involved. In [24], [25], Geroliminis and 
Daganzo studied what happen experimentally on real towns. 
• The kinetic model (Prigogine-Herman [44]) gives the evolution of the density of particles 
p(t, x, v) as a function of t, x and v, where v is the speed of particles. The model is given 
by: 

d t p + vd x p = C(p, p) , 

where C(p,p) is an interacting term that is, in general, quadratic in p and, as such, models 
the driver behavior in a simple way. From the distribution p, we can derive all the useful 
quantities, such as the average speed v(t,x) = J vp(t,x,v)dv. 
The integration of the second model is more time consuming and, therefore, not used in practice. 
The main interest of the Prigogine-Herman equation is that we can derive macroscopic laws like 
the fundamental diagram from its solution. 

The first model (the Lighthill-Whitham-Richards Model [34]) assumes the knowledge of the 
fundamental traffic diagram /. This function has been studied not only experimentally but also 
theoretically using simple microscopic models such as exclusion processes [16], [7], cellular 
automata [3], [9], or simulation of individual car dynamics. Here, we first recall a way to derive 
an approximation of this diagram on a circular road, then we present an extension to the two 
dimensional traffic (crossing roads). This derivation consists of computing the eigenvalue of a 
simple minplus linear system counting the number of vehicles entering in a road section. 

The main result obtained is the generalization of this fundamental law to the 2D case where 
the roads cross each other. In statistical physics, a lot of numerical work ([40], [12], [39], [5], [8]) 
and good surveys [9], [26] have been done on idealized towns. These works analyze numerical 
experimentations based on various stochastic models with or without turning possibilities and 
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show the existence of a threshold of the density at which the system blocks suddenly. The 
particular case of one junction is studied in [28], [29] where precise results are given for the 
stochastic case without turning possibilities. Let us mention also the attempt to derive a 2D 
fundamental diagram by Helbing [27] using queuing theory. In the saturated and unsaturated 
traffic case he derives a formula for an intersection that he extends to an area of a town that fits 
well to experimental datas. 

Here we present a deterministic model with turning possibilities, based on Petri nets and 
minplus algebra. The minplus linear model on a unique road can be described in terms of event 
graph which is a subclass of Petri nets. The presence of conflicts at junctions prevents the 
extension of this model to the 2D case. 

We propose a way to solve the difficulty by extending the class of weights used in Petri nets 
allowing negative weights. Due to such weights, the firing of a transition can consume tokens 
downstream, and the modeling power of the Petri nets is improved significantly. This possibility 
is used to model the authorization to enter in the junctions. 

The dynamics of general Petri nets allowing negative weights can be written easily but is 
neither linear in minplus algebra nor monotone. Nevertheless, in the traffic applications given 
here, dynamics (x k+1 = f(x k )) are always homogeneous of degree one (/(A ® x) = A ® f(x), 
where ® denotes the minplus multiplication that is the standard addition). For such systems, the 
eigenvalue problem (computation of A such that A <g> x = f(x)) can be reduced to a fixed point 
problem. The existence of the growth rate x — lim^a^/A; is due to the existence of a Birkhoff 
average. The quantities A and x coincide when / is also monotone, but this is not true in the 
general case. The monotone case has been studied carefully in [32], [38]. In traffic examples of 
roads with junctions, the dynamics is not monotone. 

In all the traffic examples, we study systems of roads on a torus without entries, in such a way 
that the number of vehicles remains constant in the system. Thus, we represent an idealization of 
constant densities for a more realistic system. To maintain a constant density in an open system, 
a new vehicle has to enter each time another one leaves the system. This is mathematically 
equivalent to consider circular roads. If we want to maintain a constant local density in a large 
system, it is the same problem; each time a vehicle leaves the local subsystem, another one has 
to enter. 
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We study indetail the particular case of two circular roads with one junction. For this system, a 
result on the existence of the growth rate is obtained using the nondecreasing trajectory property 
of the states starting from zero. From this result, we show that the distances between the states 
stay bounded. The eigenvalue problem can be solved explicitly. On simulation, we see that the 
eigenvalue and the growth rate generally do not coincide. However, these quantities are very 
close for any fixed density. Therefore, the simple formulas obtained for the eigenvalue give a 
good approximation of the traffic fundamental diagram. 

This definition of the fundamental diagram must not be confused with the one used by 
Daganzo [13], which is the Hamiltonian of the traffic dynamics interpreted as a dynamic pro- 
gramming equation in the simplest case of a single road and generalized in heuristic way to 
the city case. His fundamental diagram mainly coincides with our minplus dynamics, in the 
case of one street, but is observed experimentally not derived of a microscopic model as we do 
here. The dependence of the eigenvalue or the average growth rate with the density of vehicles 
in the system, which is a difficult result to obtain even in the case of two streets with one 
intersection, is not done in the Daganzo work. Moreover, for the city case, we remark also, 
on the simple example studied here, that the dynamics is not a standard dynamic programming 
equation associated to a deterministic or a stochastic control problem (since the dynamics is 
not monotone) which could be a problem, at least at the conceptual level, for the point of view 
adopted by Daganzo. 

The fundamental diagram that we obtain for a system of two roads with one junction, shows 
four phases: 

• Free phase: The density is low where the vehicles do not interact. 

• Saturation phase: The junction is saturated, but the locations downstream of the junction 
are freed when a vehicle wants to leave the junction. 

• Recession phase: The locations downstream of the junction are sometimes crowded when 
a vehicle wants to leave the junction. 

• Freeze phase: The vehicles cannot move. 

These phases are derived from a unique simple model and are not postulated to obtain different 
models subsequently analyzed. 

Preliminary results on the traffic dynamics have been presented in [18], [19], [20]. In [21], 
many developments, complementary results, and other completely solved examples can be found. 
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The theorems given here on the eigenvalue problem and the growth rate complete some of the 
main results of [21] by relaxing some assumptions and clarifying the growth rate existence 
problem. 

We show that a system theory can be developed for the class of concave polyhedral 1- 
homogeneous systems, which we construct using parallel, series, and feedback compositions. 
Moreover, the dynamics of these systems are characterized by the composition of a standard 
matrix and a minplus matrix. The sum of the entries in each line of the standard matrix is equal 
to one. We can compute easily the transformation on these two matrices corresponding to the 
three composition operations. The interest of this "system theory" can be shown by building the 
traffic dynamics of a regular town starting from three elementary systems. We do not detail the 
construction, but we show a stationary vehicle distribution obtained for a simple regular town 
with dynamics built in that way. In this case, the fundamental diagrams still present the four 
phases observed in the simple case of two roads with one junction. All the numerical simulations 
have been done using the ScicosLab software [46]. 

This paper is divided into three parts. First, we recall the basic results of minplus algebra; 
we present a system theory for polyhedral concave 1 -homogeneous systems and discuss their 
growth rate and eigenvalue problems. Second, we present Petri nets allowing negative weights. 
Third, we give applications to the computation of the fundamental traffic diagram; we consider 
the case of one circular road with minplus linear dynamics and study the case of two circular 
roads with a junction in detail; then, we explain briefly the way to build the dynamics using the 
system theory described in the first part. 

II. Minplus Algebra and extensions 
A. Review of minplus algebra 

In this section, we revisit the main definitions and properties of the minplus algebra. An 
in-depth treatment of the subject is in [4]. 

The structure K m i n = (KU {+00} , ©, <g>) is defined by the set R U {+00} endowed with 
the operations min (denoted by ©, called minplus sum) and + (denoted by ©, called minplus 
product). The element e = +00 is the zero element e © x = x, and is absorbent e © x = e. 
The element e = is the unit element e © x = x. The main difference with respect to the 
conventional algebra is the idempotency of the addition x © x = x and the fact that the addition 
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cannot be simplified a®b = c®b^a = c. This structure is called minplus algebra. We will 
call R min the completion of R min by — oo with — oo <g> e = e. 

This minplus structure on scalars induces an idempotent semiring structure on m x m square 
matrices with the element-wise minimum denoted by © and the matrix product defined by 
(A ® B) ik = min, (Aij + Bj k ), where the zero and the unit matrices are still denoted by e and 
e. We associate a precedence graph G(A) to a square matrix A where the nodes of the graph 
correspond to the columns of the matrix A and the edges of the graph correspond to the non-zero 
(t^ e) entries of the matrix. The weight of the edge going from i to j is the non-zero entry A^. 
We define the weight of a path p in Q{A), which we denote by \p\ w , as the minplus product of 
the weights of the edges composing the path (that is the standard sum of weights). The number 
of edges of a path p is denoted by \p\i. We will use the following fundamental result (see [4]). 

Theorem 1: If the graph Q(A) associated to an m x m minplus matrix A is strongly connected, 
then the matrix A admits a unique eigenvalue A G K m i n \ {e}: 



3 X e R™ in , X ^e: A®X = \®X with A 
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where C is the set of circuits of Q{A) □ 



B. A generalized matrix calculus 

In a Petri net, two kinds of operations appear. One is the accumulation of resources in 
the places, and the other is the synchronization in the transitions. The first operation can be 
modeled by addition, and the second by a min or max (a task can start at the maximum of the 
arrival instants of the resources needed by the task). Matrix notations can be generalized to such 
situations. 

We consider the set ofmxm matrices where the rows [resp. columns] are partitioned in two 
sets the standard and the minplus (here the m! first rows [resp. columns], and the last m" rows 
[resp. columns]) with entries in MU {+00, — 00}, equipped with the two operations EH and M 
defined by: 
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Since the entries are in an extension of 
multiplication table: 



we have to specify the scalar addition and 



x ±00 = ±00 x = 0, +00 (g> (— 00) = +00 — 00 = +00 . 

These choices have been made to preserve the absorption properties for the multiplication of the 
null elements of the standard algebra (0) and of the minplus algebra (e = +00). This absorption 
property is useful to model the absence of an arc in the precedence graph Q(A) associated to a 
square matrix A (defined in the same way as in the pure minplus case). 

'o 



The addition EE is associative, commutative and has the null element 



e e 



still denoted e. 



The multiplication E3 has no identity element. It is neither associative nor commutative nor 
distributive with respect to the addition. 

The main interest of this operation appears when the left operand is a vector Y = A M X, 
where X 6 i™ n is a vector and A is a m x m matrix with entries in IR m i n . the vector Y, seen 



m 
min 



as a function of X, is a set of m! standard linear forms and of m" minplus linear forms on 
with m = m! + m" . The operation Z = A M (A M X) corresponds to the compositions of these 
linear forms but these compositions do not define anymore a set of standard and minplus linear 
forms and Z/(ABA)8I. 

Moreover, when the left operand is a vector, this multiplication can be represented by the graph 
Q(A) where we have two kinds of nodes (those corresponding to the standard + operation, and 
those corresponding to © operation) and two kinds of edges (those which operate multiplicatively 
(x) and those which operate additively (®)). 

a b 



Example 1: Let us consider the graph Q(A) associated to the matrix 

one node associated to the standard algebra and one node to the minplus a 
where y and x are two vectors with two entries, means: 



(Figure 

c d 
gebra. Then y 



1 ) with 

= A^X, 



2/1 = axi + bx 2 , V2 = min(c + x 1 ,d+ x 2 ) . 



We adopt the following conventions to solve some ambiguities in the formula notations: 
- As soon as a minplus symbol appears in a set of formulas, all the operations must be 
understood in the minplus sense with the exception of the exponent that must be understood in 
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Fig. 1. The incidence graph associated to the matrix A with one ©node and one +node. 

the standard sense. For example, x a l h © 1 or \fxf- © 1 must be understood as min((a/6)x, 1) and 
not as min((a — b)x, 1). In minplus sense, a/b = a — b since it is the solution^ of b © x = a. 
The equation b © x = a means b + x = a. 

- In a minplus formula, the rational numbers (written with figures) are denoted in the standard 
algebra. For example, \x®l (instead of \flx®l) means min(0.5 + x, 1) and not min(— 1+x, 1), 
but (a/b)x © 1 means min(a — b + x, 1). 

- A non-zero element in the minplus sense means a finite element in the usual sense. Positive 
element will be used always in the standard sense. 

With these conventions the formulas are concise and not ambiguous. 

C. A generalized system theory 

We can develop a generalized system theory based on the two operations EB and IE1. We partition 
the states [resp. inputs, outputs] into two classes, the standard states [resp. inputs, outputs] and 
the minplus states [resp. inputs, outputs]. Then the dynamics can be written: 



X k+1 
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X k 
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'For readers familiar residuation theory, the minplus division used here means the standard minus operator with the convention 
previously given for infinite elements. This choice is incompatible with the residuation which chooses the smallest solution of 

b ® x — a. 
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These dynamics, denoted by S and defined by the matrices (A, B, C), associate the output signals 
(Y k ) keN to the input signals (U k ) keN - We write Y = S(U) . 

Let us note that E3 corresponds to the definition given previously when there is a permutation of 
rows and columns. This is because [X k U k ]' [resp. [X k Y k ]')] is not written in the canonical form 
since all the standard entries are not followed by all the minplus entries but by the standard states, 
the minplus states, the standard inputs [resp. outputs], and the minplus inputs [resp. outputs]. 
See the Petri net application Q for developed formulas. 

On these systems, we define the following operations: 

• Parallel Composition. Given two systems Si and S2 (with the same numbers of inputs 
and outputs), we define the system S = Si EH S2 obtained by using the same entries and 
adding^] the outputs. The dynamics of S is defined by: 



.4 



Series Composition. Given two systems Si and S2 (where the output numbers of S2 equal 
the input numbers of Si), we define the system S — Si M S2 obtained by composition of 
the two systems: S(U) = Si(S 2 (U)). Using the new state [X 1 X 2 Y 2 )', the dynamics of S 
is defined by: 
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• Feedback. Given a system S, we define the feedback system S B as the solution in Y of 
Y = S(U EH Y). Using the new state [X Y]', the dynamics of S B are defined by: 
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D. Additively homogeneous autonomous dynamic systems 

Let us now discuss one-homogeneous minplus dynamical systems which are a large class 
applied frequently when there is a conservation such as probability mass or tokens in Petri nets. 



The standard output nodes are added in the standard algebra and the minplus output nodes are added in the minplus algebra. 
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Because we accept negative entries, we obtain a generalization of two cases. One case is a 
measure not necessarily positive with a total mass equal to one. The other case is a conservation 
of tokens in a situation where negative entries may appear. Let us start with an academic example. 
The traffic modeling will give more concrete examples in the following sections. 

Example 2: Let us go back to Example [T] and assuming that a + b = 1. Adding A to each 
component of x (that we can write \<g>x) implies that the two components of y are augmented 
by A. We have 

A M (A ® x) = A <g) (A R x) . 

We say that the system is additively homogeneous of degree 1 or, more simply, homogeneous. 
Indeed, using the minplus notation, y = A M x can be written: 

Vi = (a;i) 0a <8> (z 2 )® 6 , y 2 = c®x 1 ®d®x 2 , 

which is clearly homogeneous of degree 1 in the minplus algebra as long as a+b = 1. Moreover, 
if a and b are nonnegative, then the transformation is nondecr easing. 

To simplify the notations, we will write the transformation in the following way: 

Hi = (x 1 ) a (x 2 ) b , y 2 = cxi © dx 2 . 
More generally, we say that a function / : IR™ in i-> E min is homogeneous if 

f{\®x) = A® f(x) . 

E. Eigenvalues of homogeneous systems 

The eigenvalue problem of a function / can be formulated as finding non-zero x G IR™ in and 
A G Mmin such that: 

A <g> x = f(x) . 

Since / is homogeneous, we can suppose without loss of generality that if x exists then x\ ^ e. 
The eigenvalue problem becomes: 

A = ftix/xx) , 
X2/X1 = {f2/h){x/x 1 ) , 

< 

Xn/Xi = {fn/ h){x/x\) , 
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where the division is in minplus sense, that is the subtraction. Denoting y = (x 2 /xi, • • • , x n jx\) 
and gi-i(y) = (/i//i)(0, y), the problem is reduced to the computation of the fixed point problem 
V — 9{y)- This fixed point gives the normalized eigenvector from which the eigenvalue is deduced 
by: A = fi(0,y). We note that g is a non-homogeneous minplus function of y. 

The fixed point problem does not always have a solution, but nevertheless, there are cases 
where we are able to find one: 

• / is affine in standard algebra. In this case f(x) = Ax + b. The homogeneity^] implies 
that the kernel of A — Id is not empty. When this kernel has one dimension, the eigenvalue 
is equal to A = pb, where p is the normalized (pi = 1) left standard eigenvector of A 
associated to the (standard) eigenvalue 1. But even in this case, all the standard eigenvalues 
do not have a module necessarily smaller than one, and the dynamical system may be 
unstable. We note that when all the entries of the matrix A are nonnegative, / is monotone 
nondecreasing, but when there are positive entries and negative entries, the system is not 
monotone. 

• / is minplus linear, fix) = A x. In this case the system is monotone. 

• / corresponds to stochastic control. In this case f(x) = D (Hx) where if is a standard 
matrix with rows that define discrete probability laws. Such a matrix is called a stochastic 
matrix. Then the dynamics x k+1 = f(x k ) has the interpretation of a dynamic programming 
equation associated to a stochastic control problem, and the eigenvalue is the optimal 
average cost of the corresponding stochastic control problem. Note that in this case, x k 
are components of all the dynamics which can be written 

x k+l = A^x k , with A = 

• f corresponds to stochastic games. In this case / = D\ (D 2 (Hx)) where denotes 
the maxplus product (obtained by replacing min by max in the minplus matrix product 0), 
and H is a stochastic matrix. This case corresponds to dynamic programming equations 
associated to stochastic games. In this case, / is monotone. 

3 Al — I with 1 the vector with all its entries equal to 1. 



H 

D e 
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• / has a particular triangular structure for example: 

x k+l = A®X k , 
yk+1 _ B(^ x k^ (g, yk ^ 

with B(x) additively homogeneous of degree 0. The system is not necessarily monotone. 

For such systems, it is easy to find the eigenvalue and eigenvector by applying the minplus 

algebra results. See [21], [23] for discussions and generalizations. 
In general, it is possible to compute the fixed point using the Newton's method^} This method 
corresponds to the policy iteration when the dynamic programming interpretation holds true. For 
stochastic control problems, the policy iteration is globally stable. In the game case, the policy 
iteration is only locally stable. 

We may have unstable fixed points which are not accessible by integrating the dynamics. In 
this case, the eigenvalue, computable by the Newton's method, gives no information on the time 
asymptotic of the system. When all the fixed points are unstable, we may have a linear growth 
of the state trajectories. This point is illustrated by the chaotic tent dynamics example given in 
the next section. 

F. Growth rate of homogeneous systems 

We define the growth rate, of a dynamic system x k+1 = f(x k ), where x G E^ in , by 

the common limit \im k x k /k of all the components i when this limit exists. In [32] it has been 
proved, with a special definition of connexity (satisfied for a system defined by f(x) = AM x 
as long as the graph Q{A) is strongly connected^]), that the growth rate and the eigenvalue of 
a homogeneous and monotone system exist and are equal. Let us show that chaos may appear 
and that the eigenvalue and growth rate which exist are not equal on a system which is only 
homogeneous. 

Let us consider the homogeneous dynamic system where k is the time index: 

^k+l _ k 

= (x k 2 y/{x\y®2(x k 1 y/x k 2 . 

4 we have to solve a piecewise linear system of equations 

5 Here there is an edge from i towards j in G(A) if Aj t 7^ e if it is a minplus edge or Ay, 7^ if it is a standard one. 
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The corresponding eigenvalue problem is 

{Xxi = x 2 , 
\x 2 = x\jx\ © 1x\jx 2 , 
where the minplus power exponent must not be confused with a time index. 




Fig. 2. Cycles of the tent transformation. The abscissa xm of M is a fixed point of /: xm = /(ia/)- The pair (x a , x c ) is a 
cycle of / composed of two fixed points of / o /: x c — f(x a ), x a — f(x c ) = f(f(x a ))- The triplet (a;i, X3, X5) is a circuit 
of / composed of three fixed points of / o / o /: xz = f(x%), X5 = f(xs) = f(f(xi)), Xi = f(xs) — f(f(f{xi))). 

The solution is A = y with y = x 2 jx\ satisfying the equation 

y = y 2 © 2/y 2 , 

which has the solutions y = and y = | = 0.66.... These two solutions are unstable fixed points 
of the transformation f(y) = y 2 © 2/y 2 . However, the system y k+l = f(y k ) is chaotic since / 
is the tent transform (see [6] for a clear discussion of this dynamics). In Figure-|2j we show the 
graph of x >-)■ f(x), x >-)■ f(f(x)) , x i-> f(f(f(x))), their fixed points, and periodic trajectories. 

In Figure-|3} we show a trajectory for an initial condition chosen randomly with the uniform law 
on the set {(i — 1)/10 5 , £ = 1, • • • , 10 5 }. The diagonal line in the picture is a decreasing order 
applied to the set {y k , k = 1, • • • , 10 5 }. This line shows that the invariant empirical density 
is approximately uniform. In fact, with these initial conditions, the trajectories are periodic 
with a possible long period. It has been proved that the tent iteration has a unique invariant 
measure absolutely continuous with respect to the Lebesgue measure (the uniform law on [0, 1]). 
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Fig. 3. A tent iteration trajectory (10 5 y )k=i,.,30o) an d its arrangement in decreasing order (almost diagonal line) . 



Therefore, the system is ergodic. The growth rate 



1 N 

(af-*?)/tf=^£(a?-*r: 

k=l 



can be computed by averaging, with respect to the uniform law, an increase in one step: fx (x) —x\ 
with the standard notations (that is f\(x)/xi = y with the minplus notations). Therefore x(f) — 
Jq ydy = 0.5 , for almost all initial conditions, which is different from the eigenvalues (0 and 



10 



3 J 

More generally, for a homogeneous system, we can write the system dynamics: 



x 



i fl /^ = /iWM = %*). y k+1 = g(y k 



with y\_ x = x\jx\ and = f%/ f\ for i — 2, • • • , n. As long as the y k belong to a bounded 
closed (compact) set for all k, we remark (after Kryloff and Bogoliuboff) that the set of measures: 

(where S a denotes the Dirac mass on a) is tight. Therefore, we can extract convergent subse- 
quences which converge toward invariant measures Q y o that we will call Kryloff-Bogoliuboff 
invariant measure. Then we can apply the ergodic theorem at the sequence (y k )keN- Application 
of this theorem shows that, for almost all new initial conditions chosen randomly according to 
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the Q y o, we have : 

X (f) = lim l(af - x?) = lim 1 (f^ Kv")) = J Ky)dQ y0 {y) . (2) 

It may happen that the initial condition y° is transient; therefore, it is in the attractive basin 
of Q y o not in the support of Q y o. It would be very useful to prove that y° is generic (in the 
sense of Furstenberg[31]), that is the limit exists for y°. A priori homogeneous systems do not 
have the uniform continuity property required to prove the convergence of the Cesaro means 
for all y°. The following classic example show the case of an ergodic system where with non 
generic points / : x E T 1 -> 2x E T 1 with: x° = 0.100111100000000 • • • where T 1 the torus of 
dimension 1. In Figure-|4]we see that the Cesaro does not converge. 

0.70-r 
0.65- 
0.60- 
0.55- 
0.50- 
0.45- 
0.40- 
0.35- 

0.30-| 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

200 400 600 800 1000 1200 1400 

Fig. 4. Plot of S(n) with: S(n) = A ^ ''■ 

In the case where the compact set is finite, we can apply the ergodicity results on Markov 
chains with a finite state number to show the convergence of towards Q y o. Instead of the 
subsequence convergence, this convergence proves the genericity of y°. 

This discussion is summarized by the following result. 

Theorem 2: For any additively 1-homogeneous dynamical systems x k+1 = f(x k ), x° such 
that x k jx\ stays bounded for all j and k, there exists a measure on the initial condition x° such 
that the growth rate exists for almost all the initial conditions. 
We can also see [2] for construction of invariant measures of stochastic recursions. 
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Coming back to the tent example, according to the initial value y°, the tent iterations y k stay 
in circuits or follow trajectories without circuit (possibly dense in [0, 1]). For example, assuming 
that the initial condition is y — |, the trajectory is periodic of period 2. The invariant measure 
is Q y o = + 5§_) . The growth rate is |, which is again different from the eigenvalues 
and |. Moreover, it can be shown that for all initial conditions with a finite binary development 
(this set contains all the float numbers of computers), the trajectory stays in the unstable fixed 
point after a finite number of steps. That is, for a dense set of initial conditions the invariant 
measure is <5 and the growth rate is 0. 

III. Petri Net Dynamics 

A. Autonomous Petri nets 

We present, in min-plus-times algebra, timed continuous Petri nets with weights. The weights 
can be negative and the numbers of tokens are not necessary integer (in continuous Petri nets, 
what we call tokens are in fact fluid amounts). 

A Petri net J\f is a graph with two sets of nodes (the transitions Q (with \Q\ elements) and 
the places V (with \V\ elements)) and two sorts of edges, (the synchronization edges (from a 
place to a transition) and the production edges (from a transition to a place)). 

A minplus \Q\ x \V\ matrix D, called synchronization^ matrix is associated to the synchro- 
nization edges. D qp = a p if there exists an edge from the place p E V to the transition q E Q, 
and D qp = e elsewhere, where a p is the initial marking of the place p, which is, graphically, 
the number of tokens in p. We suppose here that the sojourn time in all the places is one unit 
of time(!l 

A standard algebra \V\ x \Q\ matrix H, called production^ matrix is associated with the 
production edges. It is defined by H pq = m pq if there exists an edge from q to p, and 
elsewhere, where m pq is the multiplicity of the edge ^ 

6 Decision matrix in stochastic control. 

7 When different integer sojourn times are considered, an equivalent Petri net with a unique sojourn time can be obtained by 
adding places and transitions and solving the implicit relations. 
8 Hazard matrix in stochastic control 

'Here the multiplicity appears only with the output transition edges. The multiplicity of input transition edges is supposed to 
be always equal to one. Looking at the more general case [10], we see that we do not lose generality by doing so (the dynamics 
class obtained is the same) 
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Therefore, a Petri net is characterized by the quadruple: 



(V, Q, H, D) . 



A Petri net is a dynamic system in which the token (fluid) evolution is partially defined by the 
transition firings, saying that a transition can fire as soon as all its upstream places contain a 
positive quantity of tokens (fluid) having stayed at least one unit of time. When a transition fires, 
it consumes a quantity of tokens (fluid) equal to the minimum of all the available quantities being 
in the upstream places. Cumulating the firings done up to present time defines the cumulated 
transition firing of the transition. The firing produces a quantity of tokens (fluid) in each 
downstream place equal to the firing of the transition multiplied by the multiplicity of the 
corresponding production edge. If the multiplicity of a production edge, going from q to p, is 
negative, the firing of q consumes tokens (fluid) of p. 

A general Petri Net defines constraints on the transition firing. Denoting by q k the cumulated 
firings of transitions q E Q up to instant k, they satisfy the constraints: 



where p E V is a place of the Petri Net; q m C V [resp. q out C V] denotes set of places upwards 
[resp. downwards] the transition q; and p m C Q [resp. p out C Q] denotes the set of transitions 
upwards [resp. downwards] the place p. 

Indeed, being at time k, we know (from the firing definition of transitions) that after the firing 
(which is instantaneous), there is at least one place upstream of any transition in which no token 
entered before time k — 1. For each transition q, the equation p]) computes the number of tokens 
which have stayed at least one unit of time in each place p E q m . The equation expresses that 
at least one place is empty and the others have nonnegative numbers of tokens. 

As long as there is more than one edge leaving a place, the trajectory of the system is not 
well-defined because we do not know the path of a token leaving this place. 

In the case of a deterministic Petri net (generally called conflict free Petri net) where all the 
places have only one downstream edge, the dynamics are well-defined, meaning there is no token 
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consumption conflicts between the transitions downstream of each placej^j Then, denoting by 
Q = (Q k )qeQ,keN me vector of sequences of cumulated firing quantities of transitions, and by 
P = (p k ) p <=-p,keN the vector of sequences of cumulated token quantities arrived in the places at 
time k, we have: 



pk+l 




H 




pk+l 


def 




Q k+1 




D e 




Q k 




D <g> P k+1 



(4) 



In Equation P k+1 counts the number of tokens available (for firing in the downstream 
transitions) at time k + 1 coming from the upstream transition firings. It is obtained by summing 
the weighted firing numbers of the transitions upstream the places since the tokens are supposed 
to stay at least one unit of time in the places. 

The part about Q k+1 in Equation Q tells that the transition firing numbers at time k + 1 
are equal to the minimal token numbers in the places upstream the transitions. To obtain these 
quantities we have only to add the token numbers in the place at initial time to the numbers 
entered by the firings that are the entries of P k+1 . 

From these dynamics, we deduce the dynamics of the cumulated firing quantities by eliminat- 
ing the place variables. We deduce the dynamics of the cumulated token quantities by eliminating 
the transition variables. 

Q k+l = D® {HQ k ), P k+1 = H(D® P k ) . 

In the case of event graphs (particular deterministic Petri nets where all the multiplicities m pq 
are equal to 1 and all the places have exactly one edge upstream), the dynamics are linear in 
the minplus sense : 

Q k + l = A®Q k , 
where A q > q = a p with p the unique place between q and q'. 

B. Deterministic Petri nets 

By using negative weights and/or fixing a routing policy, it is possible to transform a Petri 
net with conflicts to a deterministic Petri net. Let us discuss these points more precisely on the 
simple system given in the first picture of Figure-[5] 

l0 In the non-deterministic case, we have to specify the rules which resolve the conflicts by, for example, giving priorities to 
the consuming transitions or by imposing ratios to be followed. As long as these rules are added, the initial non-deterministic 
Petri net becomes a deterministic one. 
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Undefined Dynamics Given Routing West p r i ority 

Fig. 5. The nondeterministic Petri net, given in the left figure, is made clear by: - choosing a routing policy: 1/2 (to read 
in standard algebra) towards q$, 1/2 towards q± in the central figure, - giving top priority to q$ against q± in the right figure 
(where the time shift given by the sticks are not anymore only in the places but also on the edges). 



The incomplete dynamics of this system can be written in minplus algebrq | 

qlql = aq\- x q\- x . (5) 

Clearly q 3 and q^ are not defined uniquely. In the following two examples, we complete the 
dynamics in two different ways. These two ways are equally useful for traffic applications as 
we will see later. 

• By specifying the routing policy (for example we choose arbitrarily that half of the total 
tokens available are given to q 3 and half to q^, see [10] for results on the general routing 
case) 

q\ = yjq*-^- 1 , <?s = • 

The minplus product of the two equations gives the constraint ([5]). 

• By choosing a priority rule (top priority to q 3 against qj^ 

The last equation implies that the initial constraint ([5]) is satisfied. We see that the negative 
weights on q\~ x and on q* are essential to express this priority. 
In the two cases, we obtain a degree one homogeneous minplus system. 

"Which means in standard algebra: 94+53 = a + q\~ 1 + q 2 ~ 1 . 

12 Which means in standard algebra: q\ = ~ (<7i _1 + <?2 _1 ) 1 Is = a + 14 ■ 

13 Which means in standard algebra: g| = a + q^ 1 + q\~ x — q^ 1 , q\ — a + gj*' _1 + q^ 1 — g| . 
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C. Input-Output Petri nets 

We can define a Petri net with inputs and outputs in the following way. We partition the 
transition set in three parts (V, Q, Z): the input set V, the state set Q and the output set Z. We 
do the same for the place set and we get the three parts (U, V, y). The inputs are the transitions 
[resp. places] without upstream edges. The outputs are the ones without output edges. Then the 
dynamics can be rewritten: 



pk+l 







A 





B 




pk+l 






AQ k + BV k 


Q k+1 




c 


£ 


D 


e 


m 


Q k 






5 pk+i Tjk+i 









E 










jjk+i 






EQ k 


Z k+l 




F 


e 


£ 


e 




yk 






F ® P k+1 



(6) 



These dynamics, denoted by S and defined by the matrices (A,B,C,D,E,F), associate the 
output signals (Y k , Z k ) km to the input signals (U k , V k ) k&N . We write (Y, Z) = S{U, V) . These 
systems can be considered as a special case of an extension of the generalized system described 
in the minplus section. In this special case, some blocks are null and an implicit part can appear. 

The three operations (parallel, series, and feedback compositions) described previously can 
be extended easily to this case. 

IV. Traffic Application 

A. Traffic on a circular road 

Let us recall the simplest model to derive the fundamental traffic diagram on a single road. 
The best way to obtain this diagram is to study the stationary regime on a circular road with 
a given number of vehicles^} We present two ways to obtain this diagram: one is by logical 
deduction from an exclusion process point of view (it shows clearly the presence of two distinct 
phases), and the other by computing the eigenvalue of a minplus system derived from a simple 
Petri net modeling of the road (this way will be extended to the case of roads with junctions). 
In the following, the road is cut in m sections which can contain at most one vehicle. 

14 We consider that the stationary regime of the circular road is reached locally on a standard road when its density is constant 
at the considered zone. 
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Fig. 6. On the top-left side, a circular road cut in sections. On the top-right side, its exclusion process, where 1 means that the 
section is occupied. On the middle, its Petri net representation where the ticks (1 time delay) in each place are not represented. 



B. Exclusion process modeling 

Following [7], we can consider the dynamical system defined by the rule 10 — > 01 applied 
to a binary word w. The word describes the vehicle positions at instant fc on a road cut 
in sections (each bit representing a section, 1 meaning occupied and meaning free, see II in 
Figure -[6]). Let us take an example: 

W! = 1101001001, w 2 = 1010100101, 
w 3 = 0101010011, w 4 = 1010101010, 
w 5 = 0101010101. 



Let us define: - the density p by the number of vehicles n divided by the number of sections 
m: p = n/m, - the flow ip(t) at time t by the number of vehicles going one step forward at time 
t divided by the number of sections. Then the fundamental traffic diagram gives the relation 
between (p(t) and p. 
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If P < 1/2, after a transient period, all the vehicle groups split off, and then all the vehicles 
can move forward without other vehicles in the way, and we have: 



If p > 1/2, the free place groups split off after a finite time and move backward without other 
free places in the way. Then m — n vehicles move forward and we have: 



Proof: This result has been proved in [7]. Let us give the idea of the proof. We only have 
to prove that after a finite number of steps, the separations of vehicles or holes appear. 

Assuming that the density is not greater than 1/2, let us look what happens at a cluster of at 
least two vehicles denoted by A. There are two cases: 

• The cluster behind A is separated from A by one hole, we have the configuration 10A0. 
Then at the next step we have 0A01. The cluster A has gone backward by one place. 

• The cluster behind is separated from A by more than one hole, we have the configuration 
00A0. Then at the next step, the size of A has been reduced by one. 

Then after a finite number of steps, bounded by the number of places, the size of each cluster 
of vehicles is reduced strictly. Indeed, the density being not greater than 1/2, and individual 
vehicles going forward, there is cluster of vehicles that must meet a cluster of at least two holes. 

For a density larger than 1/2, we follow the cluster of holes instead of vehicles. We can show 
by the same arguments that their sizes decrease or they go forward. Therefore after a finite 
number of steps the holes are separated by vehicles. ■ 

C. Event Graph modeling 

Consider the Petri net given in III of Figure-[6] which describes the same dynamics in a 
different way. In fact, this Petri net is an event graph. Therefore, the dynamics are linear in 
minplus algebra. The number of vehicles entered in the section s before time k is denoted g*. 



(p(t) = <p = n/m = p . 



(p(t) = ip = (m — n)/m = I — p . 



Theorem 3: 
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1/2 1 



Fig. 7. The fundamental traffic diagram showing the dependence of the average car flow on the car density. 



The initial vehicle position is given by the booleans a s which take the value 1 when the cell 
contains a vehicle and otherwise. 

We use the notation a — I — a. Then, the dynamics are given by: 

q k s +l = min{a_i + q k _ x , a s + q k s+1 } , 

which can be written linearly in minplus algebra: 

q k s +l = a s -iq k s _i © a s q k s+1 , 

where the index addition is done modulo m. 

Theorem 4: The average transition speed (car flow) <p depends on the car density p according 
to the law: 

if = min(p, 1 — p) □ 

Proof: This event graph has three kinds of elementary circuits: the outside circuit with 
average weight n/m; the inside circuit with average weight (m — n)/m; and the circuits on 
which some steps are made forward and then back with the average mean 1/2. Therefore, using 
Theorem [Tj its eigenvalue is 

(p = min(n/m, (m — n)/m, 1/2) = min(p, 1 — p) , 

which gives the average speed as a function of the vehicle density since the minplus eigenvalue 
is equal to \im t x,i(t)/t = <p for all i. ■ 
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D. Traffic on two roads with one junction 

Before discussing the dynamics of the town, let us study in detail the case of two circular 
roads with a junction (see the top-right side of Figure -[8]). 

A first trial is to consider the Petri net given in the middle of Figure-[8] This Petri net is not 
an event graph. It is a general non deterministic Petri Net. 




Fig. 8. A junction with two circular roads cut in sections (top-right), its Petri net simplified modeling (middle) and the precise 
modeling of the junction (top left). The ticks representing the time delays present in each place are not represented. 



We can write the dynamics of this Petri net using Equation-[3j but these equations do not 
uniquely determine the trajectories of the system. We have two places with two outgoing edges, 
a n and a n . At place a n , we have to specify the routing policy giving the proportion of cars going 
towards y 2 and the proportion going towards y 3 . At place a n , we may follow the first arrived 
the first served rule with the right priority when two vehicles want to enter in the junction 
simultaneously. Using Petri net with negative weights, we obtain the Petri Net of Figure-[8] 
where the junction is described precisely in the top-left part of the figure. In the place a n [resp 
a n+m ] are accumulated the cars going towards West [resp. towards South]. In the place a n [resp. 
d n+rn ] are accumulated the authorizations to enter in the junction from North [resp. from East]. 
These authorizations are obtained by subtracting the authorizations to enter from East [resp. 



from North] from the total number of authorizations to enter in the junction. See Section III-B 
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to have a precise description of the management of the priority rule. We remark that we have 
obtained a deterministic (conflict free) Petri Net with some negative weights. 

The dynamics describing the evolution of the vehicle number entered in section s before time 
k denoted can be obtained immediately from the Petri net of the top left corner of Figure-[8] 



using the minplus notation discussed in Section- II-B 



n k+l 



ai-iqfli © <HQi+v i ^ ~L,n,n + l,n + m, 

UnQlln+l/ Qn+m © a n-l<?n-l > 



u 



fc+1 

n+m 



an+ m grc+i/C + © a 



-lQ. 



n+m—1 ' 



(7) 



QnQn+m 



a x q\ 



9n+l — Bri+mV QnQn+m © a n+lQn+2 > 

where the entries satisfy the following constraints (written with the standard notations): 

• < dj < 1 for i = 1, ■ ■ ■ , n + m. These initial markings give the presence, 1, or absence, 0, 
of a vehicle in the road sections. However, here we see vehicles as fluid and can relax this 
integer constraint. Moreover, the transition firing cuts the tokens, and if we want to see the 
systems after some firings, there is not necessarily an integer number of tokens in a place. 
Therefore, it is better to accept real numbers of tokens belonging to the [0, 1] interval; 

• di = 1 — <2j for i 7^ n, n + m they give the initial free spaces in the places; 

• a n + a n+m < 1 the maximum number of cars in the junction is 1; 

• a n = a n+m = 1 — a n — a n + m give the free place in the junction. 

We remark that the system is homogeneous of degree 1 and that it is easy to write the dynamics 
using the generalized matrix product. Indeed, all the "monomials" appearing in the right hand 
side of System [7 (for example a/ <7*<Zn+m) are linear in the standard algebra and can be computed 
by a standard matrix product from the q vector. Then it is easy to obtain the complete right 
hand side from all the appearing "monomials" by a minplus matrix product. This matrix form 
is very useful to simulate the system in ScicosLab [46] since the minplus matrix product is 
implemented in this software. We remark also that the "monomials" like q^q^+i/Q 
negative entries in the standard matrices. Therefore the dynamics is not monotonic. 

We can prove the existence of a growth rate thanks to the following two results. 



k n+ m introduce 
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Theorem 5: The trajectories of the states of the junction dynamics (j7J), {q k )kei$> starting from 
0, are nonnegative and nondecreasing for all i □ 

Proof: Computing q 1 using the fact that all the Oj and a, are nonnegative, it is clear that 
q 1 > 0. Let us prove by induction that the trajectories are nondecreasing. It is true for k = 1. We 
suppose that it is true for k, that is q k > q k ~ l ), and we prove that it is also true for k + 1, that 
is q k+1 > q k . Rewrite §7§ q k+1 = fi(q k ) for i — 1, • • • , n + m. The functions /j for « 7^ n, n + m 
are nondecreasing. Therefore, for such an i, we have: 

Q k+1 = fM k ) > fM*- 1 ) > i , 

using first the induction hypothesis and then the dynamics definition. 
Let us prove the nondecreasing property of q n . 

• If Qn +1 = a n-lQn-l We have 

q k+1 = a n ^ql_ x > a n „ iq k z\ > fn{q k ~ V ) = q k n • 

• If it 1 = a n q k q k +1 /q k +m , using the dynamics, we have q k n+m < a^q^q^l/q*. There- 
fore, 

q k n +1 > a n q k n q k q k n+ Ja n+m q k - l q k n -^ 

which gives q k+1 > q k using the induction hypothesis and the assumption d n = a n+m . 
The nondecreasing property of q n+m is proved in the same way. If q^Xln = a n+m-iQn+m-i^ 

We have qn+ln = a n+m-l<lt+m-l — a n+m-l<?n+m-l — fn+m(<l k X ) = Qn+m • ^ Qn+m = 

a n + m Q k Q k +i/Q k+ \ we have q k+1 < a n q k q k +1 /q k +m using the dynamics. Therefore, q k n X m > 
a n+m q k +m /a n , which gives the result using a n = a n+m . U 
Theorem 6: The distances between any pair of states stay bounded. 

3ci : sup \q k - q k \ < ci,Vi, j. 
k 

Moreover 

VT, 3c 2 : sup \q k+T - q k \ < c 2 T, Vi □ 

k 

Proof: This result comes from the following inequalities (written in minplus algebra) 
deduced from the dynamics and the nondecreasing property of the trajectories: 

Q k+1 < ai-xQ k -x ,i^l,n + l 
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qi +1 < a n ^q k q k n+m < a n ^/q k +"q k n+m < a n ^j b^ 1 q\ +1 q k +m q\ +1 < a n b^q k n+m , 
with b) = <g>i =j a,. 

n k+l < a k n k < n * a k a k+m < n * n khn+m-l n k+l . n k+l ^ Un+m k 

q n+ l ^ a n+m\J q n q n +m — ^n+ m y QnQn+m — a n+m\J q n °n <?n+l Qn+1 — a n+ m O n q, n ■ 

Therefore we have: 

q n S a n -ig„_i S • • • S \px) q n + m S [Ox) a n+ m ~iq n+m -i s ■ ■ ■ S [Ox ) q n 

The result follows from these inequalities which give bounds for all the the distances between 
two states and between the same state but at different times. ■ 

By using Theorem [6] and Theorem [2j we give an existence theorem of the growth rate of the 
system ([7]). The growth rate has the interpretation of the average traffic flow. 

Theorem 7: There exists an initial distribution on (g°/g°) J=2in+m the Kryloff-Bogoljuboff 
invariant measure such that the average flow 

X = Mm q k /k, Vi , 

k 

of the dynamical system ([7]) exists almost everywhere □ 

This result is not completely satisfactory. We would like to have the existence of the growth rate 
for the initial condition g 4 ° = for all i. Nevertheless, we note that the numerical approximation 
obtained by simulation of the growth rate always exists. Indeed, in this case, the number of the 
approximated states obtained by floating number approximation is finite since they belong to a 
bounded sel[^J The existence of the growth rate can be proved using the Cesaro-convergence 
of the probability to be in a state toward an invariant measure in the case of finite Markov 
chains. At this point, it is also useful to recall that the continuous Petri net model used here is 
an approximation of a more realistic discrete Petri net. 

Theorem 8: Starting from q° = 0, if an average growth rate x exists, then x < 1/4. 
Proof: From the dynamics ([7]) we have : 

n k+1 < n n k n k I n k+1 
Hn+m — u n+mHl Hn+XI Hn i 

q\ +l < a n\j q^qn+mi 

Qn+X — a n+m\J QnQn+m- 

15 We suppose here that the approximation does not destroy the fact that the states stay in a bounded set . 
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then by summing these inequalities we get : 

- q$) + {q k n +1 - q k n ) + foffi - Q k n+ i) + - gJ+J < 1- 

Hence, when — >■ +00, and taking into account Theorem-[6j we obtain 4% < 1 



is. Eigenvalue Existence of the Junction Dynamics 

Let us consider the eigenvalue problem associated to the dynamics ([7]). It is defined as finding 
A and q such that: 

Xqi = ai_ x qi_ x © a^+i, % ^ 1, n, n + 1, n + m, 

n O'nQlQn+l I 'in+m © 0>n-lQn-l > 

(8) 

Agi = a ny /q n q n+m © , 

A(?n+1 = a n+m^ ( ln ( ln+m © <2n+lQ'n+2 j 

with: < a* < 1 for z = 1, • • • , n + m, = 1 — for i ^ n,n + m, a n + a n +m < 1 and 



The eigenvalue problem can be solved explicitly. 

Theorem 9: The nonnegative eigenvalues A as a function of the density d, written in the 
standard algebra, is given by: 



d 


< d < a 


a<d<(3 


min(/3, 7) < d < max(/3, 7) 


7 < d < 1 


A 


(l-p)d 


1/4 


(r - (1 - p)d)/(2r - 1 + 2p) 






with N = n + m, p = 1/N , r = m/N and d = ( 2^™=™ a i) /(N — 1) the density of vehicles, 
a = 1/4(1 -p),p=(r+ 1/2 - p)/2(l - p) and 7 = r/(l - p) □ 

If r > 1/2 when iV is large, the positive eigenvalue is unique and has the simple approximation: 

A ~ max <^ min I d, -, )■ (.))>. 



'4'2r- 1 

Proof: The whole proof is available in [22] (where the role of m and n are inverted). Only 
a sketch of the proof is given in this section. The proof has two parts. The first part consists of 
reducing the problem to a generalized eigenvalue problem in a four dimensional space (this part 
is easy and is not given here, it is available in [22]). The second part consists of a verification 
of the generalized minplus eigenvalue system of equations since we give explicit formulas for 
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all the eigenelements (this part is given in the Appendix, and is also available in details in [22]). 
The eigenelements have been obtained by solving explicitly the homogeneous affine systems 
with five unknowns achieving the minimum in the reduced system. Numerical simulations have 
suggested the affine system achieving the minimum. Knowing it, we had only to verify the 
inequalities proving that it achieves actually the minimum. 

Since we want the eigenvalue as a function of the density, we cannot use a numerical approach. 
Instead, we have to find explicit formulas which are piecewise affine. This is the main difficulty 
of the problem. 

After verifying that A < 1/4, by elimination of g^, i ^ l,n,n + l,n + m, thanks to the 
minplus linearity of the first equation of ([8]), we obtain the closed set of equations defining q if 
i — 1, n, n + 1, n + m; see Lemma [T] in the appendix: 

q n = (a n /X)q 1 q n+1 /q n+m © (6„/A n_1 )gi , 

q n+m = (a n+m /A 2 )gig n+ i/g n © (b m / X m ~ l )q n+1 , 

(9) 

qi = (a n /X) v /q n q n+m © (b n /\ n ~ l )q n , 

<?n+l = ( a n+m/A) qnqn+m 

where b n = ^^=1 a i * s tne number of cars in the street with priority; b n = ^^=1 ®i ^ s tne 
number of free places in the street with priority; b m = (S^tI+i 1 a * ^ s me num ber of cars in the 
street without priority; and b m = ^ ^ s me number of free places in the street without 

priority. 

Observing from numerical simulations that four phases exist, it is possible to specify their 
domains analytically, and give their physical interpretations by observing the asymptotic regimes: 

• Free moving. When the density is small, < d < a, after a finite time, all the vehicles 
move freely. 

• Saturation. When a < d < (5, the junction is used at its maximal capacity without being 
bothered by downstream vehicles. 

• Recession. When /? < d < 7, the crossing is fully occupied, but vehicles sometimes cannot 
leave the crossing because the road they want to enter is crowded. When 7 < f3, three 
eigenvalues exist on the interval [7,/?]. In this case, the system is blocked. 

• Freeze. When 7 < d < 1, the road without priority is full of vehicles. No vehicle can leave 
this road while the vehicle being in the junction wants to enter. 

arXiv: January 26 th , 2010 arXiv 



31 



Using these understanding of the phases we have been able to guess where the minimum are 
achieved and then to find the solution by solving a standard linear system of 4 equations and 5 
unknowns for each phase. In the Appendix, we show that Table [I] (written in standard algebra) 
gives the eigenvalue and eigenvector formulas as a function of the density d for each phase. 





0<d<a 


a < d < 


min(/?, 7) < d < max(/3, 7) 


7 < d < 1 


A 


(1 - p)d 


l 

4 


r-(l-p)d 
2r-l + 2p 





q n 


b n — (n— 1)A 


bn — (n — 1)A 


bn — (n — 1)A 




y n +m 


(n + 1)A — 2a n — b„ 


(n + 1)A — 2a„ — 6 n 


(n + 1)A — 2a n — b„ 


















Qn+l 


<3<n + m C*n 




4A — 1 + a n+m — a„ 





TABLE I 

The eigenvalues and eigenvectors as function of the car densities. 



In Figure-|9} we show the fundamental diagram obtained by simulation (using the maxplus 
arithmetic of the ScicosLab software [46]) for a particular relative size r of the two roads and 
the eigenvalue A given in Table [j] We see clearly the four phases described in the proof of the 
eigenvalue formula. On this figure, we see also that the growth rate and the eigenvalue are very 
close to each other at least for three phases (among four phases). 




Fig. 9. The traffic fundamental diagram x(^) when r — 5/6 (continuous line) obtained by simulation and its comparison with 
the eigenvalue A(d) given in Table [i] 
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F. Regular City Modeling. 




Fig. 10. The three elementary Petri nets with which it is possible to build the Petri net of a regular set of roads on a torus by 
composition. 



We can generalize the modeling approach used in the case of one junction to derive the 
fundamental diagram of a regular town on a torus described in Figure-[TT| (left side). We build 
the model by combining three elementary systems (representing one section (T) , a junction 
input (E), and a junction output (X)) with the composition operators described in the generalized 
system theory subsection of the minplus section. The details can be found in [21]. 

The asymptotic vehicle distribution for a small town composed of two North-South, South- 
North, East- West and West-East avenues is given in Figure-fTT] (right side). The fundamental 
diagram presents a four-phase shape analogous to the case of two roads with one junction. In 
this more general case, the role of the non-priority road is played by a circuit of non-priority 
roads which blocks the whole system when the circuit is full. 

V. Conclusion 

The Petri net modeling of traffic in a regular town can be done thanks to the introduction of 
negative weights on output transition edges. The dynamics have a nice degree one homogeneous 
minplus property but are no longer monotone. This loss of monotonicity implies that the eigen- 
value and the growth rate are no longer equal. Experimental results show that they are close in 
the case of two roads with one crossing where we are able to solve explicitly the eigenvalue 
problem and compute numerically the growth rate. The fundamental traffic diagram, which gives 

arXiv: January 26 lh , 2010 arXiv 



33 




the dependence of the average car-flow (given by the growth rate) on the density, presents four 
phases that have traffic interpretations. 

This set of 1 -homogeneous minplus systems seems to be a good class of systems that we 
can describe by two matrices, one in the standard algebra and one in the minplus algebra. The 
standard compositions of these systems can be easily described in terms of these two matrices. 
These compositions of simple systems are useful to build dynamics of large systems like that 
of regular town traffic. 

In a forthcoming and more traffic-oriented paper, we will describe further the traffic interpre- 
tation of the four phases valid also for more general systems like regular towns. The influence 
on the fundamental diagram of the traffic control, using signal lights, will also be studied. 

VI. Appendix 

Lemma 1: The eigenvalue problem ([8]) of size N can be reduced to the eigenvalue problem ([9]) 
of size 4. 

Proof: First let us verify that the solutions of ^ verify A < 1/4. Indeed ^ imply that 
Mn+m < a n+m q 1 q n+1 1 (Xq n ) , Aft < a ny /q n q n+m and \q n+1 < a n+my /q n q n+m . Multiplying these 
three inequalities we obtain A 4 < a n+m a n a n+m = 1. 

To obtain the result we have to eliminate % ^ 1, n, n + l,n + m, that is to solve a linear 
minplus system which has a unique solution as soon as A < 1/2. Indeed the loops of the 
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precedence graph associated to the linear system has all its loops positive when A < 1/2. 
Moreover we can compute explicitly its solution. For i = 2, n — 1 we have : 



Qi = qi 



i-l 



WA) 



,3=1 



n-1 



(In 



for i = n + l,n + m — 1 : 









n+m— 1 






(g) (%/A) 




(g) fe/A) 


Qn+m 




J=n+1 









Using this explicit solution in the four last equations of ([8]) we obtain the reduced system 



To verify the results given Table |I} let us rewrite the system ([9]) with simplified notations : 

U = q n ,V = q n+m , X = q u Y = q n+1 , g = b n , h = b m , k = a n , I = a n+m , wn' = n - 1, and 
m' = m — 1. 

U = kXY/XV © gX/\ n ' , 

V = kXY/X 2 U © hY/X m ' , 

(10) 

X = kVUV/X © gU/X n ' , 

Y = lVUV/X © hV/X m ' , 

with k = 1 — k — I >0, g = n' — g > and h = m' — h > are the free places in the crossing 
and the two roads. 



Lemma 2: The eigenvector ([/, V, X, Y) of the minplus nonlinear system (10) is given (using 
minplus notations) by : 





< d < a 


u<d<f3 


min(/3, 7) < d < max(/3, 7) 


7 < d < 1 


u 


g/x n ' 


g/x 1 ' 


g/x n ' 


kh 


V 


\ n '+ 2 /k 2 g = hl/k\ m ' 


l\ n '- 2 /k 2 g = \ n '+ 2 /k 2 g 


X n '+ 2 /k 2 g 


e/k 2 kh 


X 


e 


e 


e 


e 


Y 


l/k 


l/k 


X 2+n '- m 'h/k 2 g = XH/lk 


e/k 2 k 



(11) 



where the eigenvalue A of the minplus nonlinear system (10) has been given in Theorem [9] □ 
Proof: Let us consider the four density regions corresponding to the four columns of this 
table : 

1) < d < a : the first column of Table [j] is solution of the standard linear system : 

U = gX/X n ', V = hY/X m ', X = kVUV/X , Y = iVUV/X , (12) 
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which is itself a solution of System dTOb since : 



a) kXY/XVU = kkl/X 3 = 1/A 3 > (since A < 1/4 indeed, once more, using (10) we 
have VXY < kklVXY/X 4 ), 

b) kXY/X 2 VU = 1/A 4 > , 

c) gU/XX n ' = (1/A 2 ) n ' > , 

d) hV/YX m ' = hh/X m ' = (1/A 2 )™' > . 

Moreover, multiplying the 4 equalities of (12) we obtain \ n '+ m '+ 2 = ghlk which gives the 



value of A given in Table [TJ 
2) a < d < (3 : the second column of Table [I] is solution of the standard linear system : 



U = gX/X n \ V = kXY/X 2 U, X = kVUV/X, Y = lVUV/X, (13) 



which is itself solution of System (10) since : 

a) kXY/XVU = 1/A 3 > , 

b) hY/VX m ' = hglk/X n ' +m ' +2 = d m ' +n ' +1 /(mn) 1 ^ > since d > a , 

c) gU/XX n ' = gg/X 2n ' = (1/A 2 ) n ' > , 

d) hV/YX m ' = m'X n ' +2 - m ' /hkgl = m'(2n' /m') 1 / 4 /d m ' +n ' +1 > since d < f3 



Moreover using the equality giving two expressions for the value of V in ( 1 1 ) we obtain 

A = 1/4. 

3) min(/3, 7) < d < max(/3, 7) : the third column of Table [j] is solution of the standard linear 
system : 

U = gX/X n ', V = cXY/X 2 U , X = kVUv/X , Y = hV/X m ' , (14) 



which is itself solution of System (10) since : 

a) kXY/XVU = A > , 

b) hY/VX m ' = hh/X 2m ' = (1/A 2 )™' > since A < 1/4 , 

c) gU/XX n ' = gg/X 2n ' = (1/A 2 ) n ' > , 

d) iVUV/XY = 1/A 4 > . 

Moreover using the equality giving two expressions for the value of Y in ( fTl) we obtain 
X 2+n '~ m ' m' = klgh = (m' + n' + l)d which gives the value of A given in Table [ij 
4) 7 < d < 1 : the fourth column of Table [I] is solution of the standard linear system : 



U = kXY/XV , V = kXY/X 2 U , X = eVUV/X , Y = hV/X m ' , (15) 
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which is itself solution of System (10) since : 

a) gX/\ n 'U = g/kh = ghkl/m = d n+m - 1 /m > (since d > 7) , 

b) hY/\ m 'V = hh = w! > , 

c) gU/XX n ' = gkh > , 

d) iVUV/XY = Ikk = 1 . 



Moreover the compatibility of first two equalities of (15) implies that A = 0. 
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